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Abstract 

A search for pair-produced scalar fermions under the assumption that i?-parity is 
not conserved has been performed using data collected with the OPAL detector at LEP. 
The data samples analysed correspond to an integrated luminosity of about 610 pb^^ 
collected at centre-of-mass energies of y/s = 189—209 GeV. An important consequence of 
i?-parity violation is that the lightest supersymmetric particle is expected to be unstable. 
Searches for i?-parity violating decays of charged sleptons, sneutrinos and squarks have 
been performed under the assumptions that the lightest supersymmetric particle decays 
promptly and that only one of the i?-parity violating couplings is dominant for each 
of the decay modes considered. Such processes would yield final states consisting of 
leptons, jets, or both, with or without missing energy. No significant signal-like excess of 
events has been observed with respect to the Standard Model expectations. Limits on the 
production cross-sections of scalar fermions in i?-parity violating scenarios are obtained. 
Constraints on the supersymmetric particle masses are also presented in an i?-parity 
violating framework analogous to the Constrained Minimal Supersymmetric Standard 
Model. 
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1 Introduction 



In Supersymmetric (SUSY) [1] models, each elementary particle is accompanied by a supersym- 
metric partner whose spin differs by half a unit. Many of the searches for these supersymmetric 
particles ( "sparticles" ) are performed within the Minimal Supersymmetric extension of the 
Standard Model (MSSM) [2], assuming it!-parity conservation. i?-parity [3] is a multiphca- 
tive quantum number defined as Rp — (_ i)25+3B+l T^i^Q^e S, B and L are the spin, baryon 
and lepton numbers of the particle, respectively. i?-parity discriminates between ordinary and 
supersymmetric particles: Rp = +1 for Standard Model particles and Rp = -1 for their super- 
symmetric partners. i?-parity conservation implies that supersymmetric particles are always 
pair-produced and always decay through cascade decays to ordinary particles and the lightest 
supersymmetric particle (LSP). In this context, the LSP is often assumed to be the lightest 
neutralino, Xii which is then expected to be stable and to escape detection due to its weakly 
interacting nature. The characteristic signature of the supersymmetric i?-parity conserving 
decays is therefore missing energy. 

If i?-parity is violated, sparticles can be pair-produced but may decay directly to Standard 
Model particles. Scalar neutrinos could also be singly-produced via an s- or t-channel sneutrino 
exchange [4, 5] . In the present paper, the possible direct manifestations of i?-parity breaking 
couplings via pair-production of sparticles and their subsequent i?-parity violating decays to 
Standard Model particles are studied. The signatures sought in the analyses of this paper 
therefore differ from the missing energy signatures of i?-parity conserving processes. 

i?-parity violating interactions of the particles of the MSSM are parametrised with a gauge- 
invariant superpotential that includes the following Yukawa coupling terms [6]: 

Wrpv — XijkLiLjEk + \jj^LiQjDk + Xj^^^UiDjDk + eiLiH2, (1) 

where i,j, k are the generation indices of the superfields L, Q, E, D and U. L and Q are lepton 
and quark left-handed doublets, respectively. E, D and U are right-handed singlet charge- 
conjugate superfields for the charged leptons, down- and up-type quarks, respectively. H2 is 
the Higgs doublet superfield with a weak hypercharge Y — 1. The last term in (1) which mixes 
the lepton and Higgs superfields is not considered in this paper. The interactions corresponding 
to these superpotential terms are assumed to respect the gauge symmetry SU(3)c x SU(2)l x 
U(1)y of the Standard Model. The Xijk are non-vanishing only if i < j, so that at least two 
different generations are coupled in the purely leptonic vertices. The A ijk are non-vanishing 
only for j < k. The A and A couplings both violate lepton number conservation and the A 
couplings violate baryon number conservation. There are nine A couplings for the triple lepton 
vertices, 27 A couplings for the lepton-quark-quark vertices and nine A couplings for the triple 
quark vertices. There arc therefore a total of 45 i?-parity violating couplings. 

In the Constrained Minimal Supersymmetric Standard Model (CMSSM), which assumes 
a common gaugino mass, mi/2, and a common sfermion mass, rriQ, at the Grand Unification 
(GUT) scale, all sparticle masses and i?-parity conserving couplings are completely determined 
by mo and a set of three parameters: M2, the SU(2) gaugino mass parameter at electroweak 
scales^, II, the mixing parameter of the two Higgs doublets and tan/3 — V2/V1, the ratio of the 
vacuum expectation values for the two Higgs doublets. 

^Mi, the U(l) gaugino mass at electroweak scales, is related to M2 by the usual gauge unification condition: 
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Prom low energy experiments, there exist several upper bounds on the i?-parity violating 
Yukawa couplings, A, A and A . A list of upper limits on individual couplings can be found in [5, 
7-15]. The experimental limits on processes such as proton decay imply that it is reasonable 
to assume that no more than one of the 45 couplings is significantly different from zero [6] . 

This paper describes searches for ii-parity violating decays of pair-produced scalar fermions 
( "sfermions" ) , such as the charged and neutral scalar leptons {i and i> respectively) , scalar top 
quark (t) and scalar hght quarks (q) . The results of similar searches performed with the OPAL 
data at centre-of-mass energies up to 183 GeV have already been pubhshed [16]. 

Two different scenarios are probed. In the first scenario, the decays of sfermions via the 
lightest neutralino, Xit ^-re considered, where Xi is treated as the LSP and assumed to decay 
via an i?-parity violating interaction. These are denoted "indirect decays". SUSY cascade 
decays via particles other than the LSP are not considered; however, when calculating limits 
in the i?-parity violating framework analogous to the CMSSM, these cascade decays are taken 
into account by assuming a zero efficiency for their detection. In the second scenario, "direct" 
decays of sparticles to Standard Model particles are investigated. In this case, the sparticle is 
assumed to be the LSP, such that i?-parity conserving decay modes do not contribute. 

In both the direct and indirect decay scenarios, it is assumed that only one of the 45 Yukawa 
couplings is significantly non-zero at a time. Moreover, only values of the A-couplings (A, A , 
A") larger than O{10-^) are relevant to this analysis. For smaller couplings, the hfetime of 
sparticles would be sufficiently long to produce a secondary decay vertex clearly detached from 
the primary vertex, or even outside the detector. These topologies have not been considered in 
this paper. Decays outside the detector have been treated in [17]. 

Results published by the other LEP Collaborations can be found in [18]. Results have also 
been obtained by the CDF and DO Collaborations [19,20] at the Tevatron and by the HI and 
ZEUS Collaborations [21,22] at HERA. 

The production and i?-parity violating decays of £, u via A and A', t via A' and q via A" 
are described in Section 2, together with the possible signal topologies resulting from these 
processes. Short descriptions of the OPAL detector and of the data samples used are presented 
in Section 3. The signal and background Monte Carlo simulations used in the different analyses 
are described in Section 4. Sections 5 to 9 describe the specific analyses optimised to search for 
i?-parity violating processes. The physics interpretation is given in Section 10 which presents 
cross-section limits and interpretations in the CMSSM. A summary is given in Section 11. 

2 Sparticle production and decays 

This section briefly describes the production and decay modes of different sfermion species. 
A detailed description of the various decay channels, decay widths and assumptions can be 
found in [16]. The decay modes studied, resulting from non-zero A, A and A couplings, are 
presented in Table 1. Table 1 also summarises the production and decay mechanisms as well 
as the coupling involved in each decay and the final state topologies. In the indirect decays, 
the particles resulting from the x° decay are shown in parentheses. The sneutrino indirect 
a' decays followed by — > £qq, leading to topologies with jets and one or more leptons, are 
not investigated because the branching ratio into these final states is not significant. In the 
available phase space and for the CMSSM parameter space relevant to these analyses, the decay 
Xi — i^qq is preferred with respect to the decay Xi ~^ ^QQ because of the light sneutrino. 



5 



Production and Decay 



Coupling 



Topology 



A direct 



2£ + E„ 



vv 



£-£+ e+r 



X indirect 
A direct 



A f + E ■ 



A indirect 



6e + E„ 



tltl 
ti!i 
tltl 



e'^q e q 
fi+q fi-q 
T^q T~q 



A direct 

a' direct 
a' direct 



2 e + 2 jets 
2 /X + 2 jets 
2 r + 2 jets 



i-{i^qq') i+{i^qq') 
£-(v qq) e+{£^qq') 
qq) qq) 



A indirect 
a' indirect 
a' indirect 



> 2£ + jets 

> 2£ + jets 

> 2£ + jets 



vv vxi i^Xi 



v{^v qq) v{^v qq) 



A indirect 



4 jets + E„ 



vv — > 
qq^ 



qq qq 
qq' qq' 
qq' qq' 



A direct 
a' direct 
a" direct 



4 jets 
4 jets 
4 jets 



Table 1: List of production and decay mechanisms of the channels that are covered by the 
various analyses described in this paper. The couplings involved and decay type together with the 
resulting topologies are given in the second and third columns, respectively. The corresponding 
section number, where details of the relevant analysis are presented, is indicated in the last 
column. 
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3 Data samples and OPAL detector 



The analyses presented in this paper use the data collected at LEP from centre-of-mass energies 
{y/s) between 189 GeV and 209 GeV, including the data obtained at the highest energies at- 
tained by LEP in the year 2000. The total data sample corresponds to an integrated luminosity 
of about 610 pb~^. The data set collected at energies > 206 GeV has a luminosity weighted 
average centre-of-mass energy of 206.6 GeV. The results presented here are then combined with 
results obtained using data at lower energies [16]. The various integrated luminosities of the 
data taken between 1998 and 2000 are summarised in Table 2. The integrated luminosity is 
measured with a precision of 0.3—0.4% [23]. These uncertainties are taken into account in the 
limit calculation in all channels considered. 



Year 

(GeV) 
jC (pb-i) 


1998 


1999 


2000 


189 
179 


192 
29 


196 
73 


200 
74 


202 
37 


< 206 
87 


> 206 
133 



Table 2: Integrated luminosities recorded with the OPAL detector between 1998 and 2000. Due 
to different requirements on the operation of the OPAL suhdetectors the precise integrated lu- 
minosity differs from one analysis to another. 

A complete description of the OPAL detector can be found in Ref. [24], and only a brief 
overview is given here. 

The central detector consisted of a system of tracking chambers providing charged particle 
tracking over the polar angle^ range of | cos^| < 0.96 inside a 0.435 T uniform magnetic field 
parallel to the beam axis. The central tracking system was composed of a large volume jet 
chamber containing a high precision drift chamber and a two-layer silicon microstrip vertex 
detector [25], and was surrounded by a set of z chambers measuring track coordinates along 
the beam direction. A lead-glass electromagnetic calorimeter located outside the magnet coil 
covered the full azimuthal range with excellent hermeticity in the polar angle range of | cos^| < 
0.82 in the barrel region and 0.81 < |cos^| < 0.984 in the endcap regions. The magnet 
return yoke was instrumented to comprise the barrel and the endcap sections of the hadron 
calorimeter. With the hadron pole tip detectors included, the hadron calorimeter covered 
the region with |cos^| < 0.99. Four layers of muon chambers covered the outside of the 
hadron calorimeter. A system of forward detectors close to the beam axis completed the 
geometrical acceptance down to 24 mrad, except for the regions where a tungsten shield was 
present to protect the detectors from synchrotron radiation. The forward detectors included 
a lead-scintillator sandwich calorimeter, a ring of lead-scintillator modules, scintillating tile 
counters [26] and the silicon-tungsten luminometer [27] located on both sides, a few meters 
away from the interaction point. 

The tracks and calorimeter clusters used in the analyses were required to satisfy minimum 
quality criteria. These were slightly different for the analysis of purely leptonic final states and 
for the analyses with jets, either with or without leptons. For leptonic final states, a track was 

^Thc OPAL coordinate system is defined so that the z axis is in the direction of the electron beam, the x 
axis points towards the centre of the LEP ring, and 9 and </> are the polar and azimuthal angles, defined relative 
to the +z and +x axes, respectively. The radial coordinate is denoted as r. 
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required to have at least 20 hits in the wire tracking chambers, the first of which was required 
to be at a radius of less than 75 cm, and to have a momentum transverse to the beam direction 
larger than 0.1 GcV. Its distance of closest approach to the interaction point in the transverse 
plane (do) had to be less than 1 cm and the distance along the beam direction less than 
40 cm. Clusters of energy deposited in the barrel region of the electromagnetic calorimeter 
were taken into account in the analyses if their energy was larger than 0.1 GeV. Similarly, 
endcap electromagnetic energy clusters were taken into account in the analyses if their energy 
was larger than 0.2 GeV and if the cluster was composed of at least two adjacent lead blocks. 
Hadron calorimeter energy clusters had to contain an energy in excess of 0.6 MeV in the barrel 
and endcap regions, and 3 GeV in the pole tips. For leptons with jets and for purely hadronic 
channels, tracks in the central detector and energy clusters in the calorimeters were required 
to satisfy the quality criteria employed in [28]. 

In calculating the total visible energies and momenta of events and individual jets, correc- 
tions were applied to prevent double-counting of energy in the case of tracks and associated 
energy clusters. For the leptonic channels, an algorithm based on a global approach described 
in [29] was used. For jets with leptons and for purely hadronic channels, an algorithm matching 
tracks and calorimeter energy clusters, described in [28], was employed. 

4 Monte Carlo simulation 

Monte Carlo samples were generated corresponding to charged slepton, sneutrino, stop and 
squark pair-production processes as well as to Standard Model processes. All generated events 
were processed through the full simulation of the OPAL detector [30], and the same analysis 
chain was applied to simulated events as to the data. 

All sparticle production and decays considered in this paper, except for stop production de- 
scribed below, were simulated at various centre-of-mass energies between 189 GcV and 209 GeV 
and for various sfermion masses between 45 and 103 GeV using the Monte Carlo program 
SUSYGEN [31]. All MC samples of signal events contain 1000 events per mass point. 

Charged and neutral slcptons decaying via a A coupling were simulated at = 189 and 
205 GeV. Charged sleptons decaying indirectly via a A coupling were simulated at = 183 
and 205 GeV. The selection efficiencies were interpolated (or extrapolated) for the centre-of- 
mass energies where no samples were generated. A few test samples were produced to test the 
validity of the interpolation. 

For charged and neutral slcptons decaying directly via a A' coupling, production was done 
at ^/s = 189, 196, 200 and 206 GeV. Squark pair-production samples were generated at ^/s = 
189, 192, 196, 200, 202 and 206 GeV. In both cases, the signal samples with the closest ^/s 
were used for the centre-of-mass energies where no samples were simulated. 

For the indirect decays, events were produced with three sparticle masses (mj) and the 
Xi mass was fixed such that Am — — m^o — mi^/2. Samples were also simulated at the 
kinematic limit with Am — mi — m^o = 5 GeV to account for changes in the event topologies. 
The values of Am were chosen to cover a wide range of event topologies for a limited number 
of Monte Carlo events. For charged sleptons decaying indirectly with a A and a A coupling, 
selection efficiencies were also obtained using samples generated with Am equal to the mass of 
the lepton plus 0.1 GeV in order to cover the small Am region. To study the selection efficiency 
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dependency on large Am, samples of sneutrinos decaying indirectly with a A and a A' coupling 

were also generated with m^o = 10, 25 and 40 GeV, respectively for rric, < 50 GeV, 70 < < 

80 GeV and 90 GeV < mj>. For charged slepton indirect decays via A', events were simulated 
for each lepton flavour corresponding to the first index of A . The quark flavours corresponding 
to the second and third indices of A were fixed to the first and second generation, with the 
exception of a few additional samples containing bottom quarks for systematic checks. All 
possible sparticle pair-production processes were simulated only at 183 GeV [16] in order to 
determine the A-like couplings which give the lowest and the highest selection efficiencies. At 
higher centre-of-mass energies, only the decays governed by these two couplings were simulated 
and considered when computing results. 

Charged and neutral sleptons decaying directly via a A' coupling were generated using the 
Lund string fragmentation scheme of PYTHIA [32,33]. In this scheme, the final state quarks 
are evolved according to the "Lund Parton Shower" model. A colour string is formed between 
colour singlet qq states. The two quarks arc subsequently evolved in time and space, emitting 
gluon radiation between the two quarks. The jets are relatively soft owing to the gluon radiation 
between the quark states. In the case of squark production, only the independent fragmentation 
model was implemented at the generator level. In this scheme, the final state quarks are treated 
as independent particles with respect to each other. There is no QCD radiation between the 
final state quarks, and the hadronization process produces relatively hard jets. For systematic 
studies additional smuon samples were generated by the independent fragmentation model. 

For a stop decaying via the A' coupling into a quark and a lepton, all nine combinations 
of quark and lepton flavours in the final state were generated at ^/s = 189, 192, 196, 200, 
202 and 207 GeV. Since stop masses below 73 GeV were excluded in [16], only stop masses 
between 65 and 103 GeV were considered here. The production and decay of the stop were 
simulated as described in [34] . The stops were hadronised to form colourless hadrons using the 
Peterson fragmentation function [35], and associated fragmentation particles according to the 
Lund string fragmentation scheme. For the decay, a colour string was stretched between the 
spectator quark and the quark from the stop decay. Further hadronisation was also done using 
the Lund scheme. The left-right stop mixing angle, 9i , was set to zero. The dependences of 
the detection efficiencies on the fragmentation function, on the mixing angle and on the Fermi 
motion were evaluated by generating additional Monte Carlo samples with the appropriate 
parameters varied. 

The main sources of background arise from Standard Model two-photon processes (c'^^c" —> 
e+e~77 e"'"e~X), and from four-fermion (c+e~ ffff) and two-fermion (lepton-pair 
e+e~ — > (Z/j)* — >• /r and multi-hadronic e+e^ — > (Z/j)* — > qq) final states. For two-photon 
processes, the PHOJET [36] and HERWIG [37] generators were used to simulate hadronic final 
states. The Vermaseren [38] generator was used to estimate the background contribution from 
all two-photon c~^c~i~^£~ final states. Four-fermion final states, other than two-photon events, 
were simulated with grc4f [39], which takes into account all interfering four-fermion diagrams 
with the exception of the multiperipheral two-photon processes. An e'^e^£^£~ four-fermion 
sample generated with grc4f v2.2 with all interfering four-fermion diagrams including the mul- 
tiperipheral two- photon processes was compared with the standard combination of grc4f v2.1 
and Vermaseren events for a systematic cross-check. For final states with hadronic jets, system- 
atic studies were performed using KORALW [40] , which internally uses grc4f matrix elements. 
For purely hadronic final states, additional special samples generated with KandY [41], a spe- 
cial version of KORALW running concurrently with YFSWW3 [42], were also used. For the 
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two-fermion final states, BHWIDE [43] was used for the ee(7) final state and KORALZ [44] 
and KK2F [45] for the jiji and the rr states. The multi-hadronic events, qq(7), were simulated 
using PYTHIA [32] and KK2F. 

Final states with six or more primary fermions were not included in the background Monte 
Carlo samples but they are expected to make only a negligible contribution to the background. 

5 Multi-lepton final states 

This section describes the searches for purely leptonic final states resulting from pair-production 
of neutral or charged slcptons, with subsequent direct or indirect A decays (sec Table 1). The 
inefficiencies and systematic uncertainties associated with each analysis are estimated in a 
similar way for all multi-lepton final state searches, and are presented in Section 5.5. The 
selections of purely leptonic final states use the data collected at centre-of-mass energies between 
189 GeV and 209 GeV corresponding to an integrated luminosity of 613.6 pb~^. 

5.1 Event and lepton selection 

Multi-hadronic, cosmic and Bhabha scattering events were vetoed [46]. A loose high multiplicity 
veto was also apphed: events were rejected if the total multiplicity of tracks in the central 
detector and energy clusters in the electromagnetic calorimeters satisfying the quality criteria 
described in Section 3 was greater than 26. At the preselection level, it was also required that 
the ratio of the number of such "good" tracks as defined in Section 3 to the total number of 
reconstructed tracks be greater than 0.2 in order to reduce backgrounds from beam-gas and 
beam-wall events. At least two good tracks were required. 

Only tracks with | cos6'| < 0.95 were considered for lepton identification. Tracks resulting 
from photon conversion were rejected using the algorithm described in [47]. A track was 
considered "isolated" if the total energy of other charged particle tracks within 10° of the 
lepton candidate was less than 2 GeV. In these searches, a very loose lepton identification is 
sufficient. A track was selected as an electron candidate if one of the following three criteria 
was satisfied: (i) the output probability of the neural net algorithm described in [48] was larger 
than 0.8; (ii) the electron selection algorithm described in [49] for the barrel region or in [50] for 
the endcap region was satisfied; (Hi) 0.5 < E/p < 2.0, where p is the momentum of the electron 
candidate and E is the energy of the electromagnetic calorimeter energy cluster associated 
with the track. The simple criterion (Hi) identifies most electrons. The algorithms (i) and 
(ii) are based predominantly on E/p, dE/dx and shower shape information, and are used to 
increase the electron selection efficiency. These algorithms are also used in the lepton plus 
jets analyses. A track was selected as a muon candidate according to the criteria employed in 
opal's analysis of Standard Model muon pairs [51], that is, the track had associated activity 
in the muon chambers or hadron calorimeters or it had a high momentum but was associated 
with only a small amount of energy deposited in the electromagnetic calorimeter. Identified 
electrons or muons were also considered as tau candidates. Taus were selected by requiring that 
there be at most three tracks within a 35° cone. The invariant mass computed using all good 
tracks and electromagnetic energy clusters within the above cone had to be less than 3 GeV. 
For muon and electron candidates, the momentum was estimated from the charged particle 
track momentum measured in the central detector, while for tau candidates, the momentum 
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was estimated from the vector sum of the measured momenta of the charged particle tracks 
within the tau cone. 

For the two- and six-lepton final states, the large background from two-photon processes was 
reduced by requiring that the measured amount of energy deposited in each silicon tungsten 
calorimeter, in each forward calorimeter, and in each side of the gamma-catcher be less than 
5 GeV ("forward detector vetoes"). 

5.2 Final states with two leptons and missing energy 

Final states with two charged leptons and missing energy may result from direct charged slepton 
decays via a A coupling. The analysis was optimised to retain good signal efficiency while 
reducing the background, mainly due to W+W" — > i~^i/£~i/ events and two-photon processes. 

The following criteria were applied in addition to those described in Section 5.1. 

(1) Events had to contain exactly two identified and oppositely-charged leptons, each with a 

transverse momentum with respect to the beam axis greater than 2 GeV. 

(2) The background from two-photon processes and "radiative return" events (e"'"e~ — >■ Z7, 

where the 7 escapes down the beam pipe) was reduced by requiring that the polar angle 
of the missing momentum, 6*111188, satisfy | cos^missl < 0.9. 

(3) To reduce further the residual background from Standard Model lepton pair events, it was 

required that m^is/ ^ < 0.80, where mvis is the invariant mass derived from the measured 
energy and momentum of all the tracks and electromagnetic energy clusters observed in 
the event, hereafter referred to as event visible mass. 

(4) The acoplanarity angle^ (0acop) between the two leptons was required to be greater than 10° 

in order to reject Standard Model leptonic events, and smaller than 175° in order to reduce 
the background due to photon conversions. The acoplanarity angle distribution is shown 
in Figure 1 after cuts (1) to (3). Figure 1 demonstrates the good discriminating power 
of this quantity. There is a disagreement between data and Standard Model expectation 
which is characteristic of the presence in the data of un-modelled two-photon events with 
low acoplanarity. The acoUinearity angle^ (6'acoi) was also required to be greater than 10° 
and smaller than 175°. 

(5) It was required that there be no electromagnetic energy cluster with an energy larger than 

25 GeV and no associated track and no hadronic energy cluster with an energy larger 
than 10 GeV and no associated track. 

(6) In order to reduce the background from two-photon events with high transverse momentum, 

complementary cuts were applied for different regions of acoplanarity. For events with 
small acoplanarity (0acop < 1-2 rad), a cut was made on a™'^^, the component of the 
missing momentum vector perpendicular to the event thrust axis in the plane transverse 

^The acoplanarity angle, ^acop, is defined as 180° minus the angle between the two lepton momentum vectors 
projected into the x — y plane. 

^The acoUinearity angle, ^acoh is defined as 180° minus the 3-dimensional space-angle between the two lepton 
momentum vectors. 
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Figure 1: Final states with two leptons and missing energy: distribution of the acoplanarity 
angle after applying cuts (1) to (3). The dashed histogram shows signal Monte Carlo events for 
direct decays ofe via A122 with mg = 78 Ge V. Data are shown as points and the sum of all Monte 
Carlo background processes is shown as a solid line. Data used for this plot were collected in the 
year 2000 at a luminosity weighted average centre- of -mass energy of 206.6 GeV. The simulated 
signal events have arbitrary normalisation. The arrows point into the region accepted by the 
cuts. The disagreement between data and Standard Model expectation, visible in a few bins, is 
due to the presence of un-modelled two-photon events in the data. 
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to the beam axis. For events with large acoplanarity (0acop > 1-2 rad), a cut was made on 
pmiss^ the missing transverse momentum. In the low acoplanarity region, a™^^^ provides a 
good discriminating power since it is less sensitive than p™'*^*^ to the presence of neutrinos 
from tau decays or poorly measured tracks. However, in the large acoplanarity region, 
Qiniss ^Qgg j^Q^ discriminate sufficiently between signal and background and pf^^^ is used. 
A further subdivision was applied according to the event visible energy, Eyis, which is 
defined as the scalar sum of the energy of all tracks and electromagnetic energy clusters 
of the event. This selection provides a good discrimination between events with genuine 
missing energy and events which have missing energy due to secondary decays such as 
T decays. It is optimized to retain good signal efficiency also in the case of very small 
visible energy: 

1. Cuts in the small acoplanarity region: 

(a) If ^vis/v^ < 0.06, then af^'/y/^ > 0.015, 

(b) If ^vis/v^ > 0.06, then af'^y/^ > 0.0175. 

2. Cuts in the large acoplanarity region: 

(a) If ^vis/v^ < 0.06, then pf''/^ > 0.025, 

(b) If ^vis/v^ > 0.06, then pf^'/y^ > 0.035. 

(c) The event transverse momentum computed without the hadron calorimeter was 
required to be larger than 0.02 Xy^. This cut rejects events with a large mea- 
sured transverse momentum due to fluctuations in the hadron calorimeter. 

(7) The background from two-photon processes and W"'"W~ production was further reduced 
by categorising the events in different classes according to the flavour of the leptons in the 
final state. The flavours of the two identifled leptons were required to match the flavours 
of the leptons expected for each signal class. Events were selected by applying cuts on 
the momentum of the two leptons and on the visible energy: 

• eueu selection: E^is/y/s < 0.8 and PmaxZ-E-beam > 0.1, where Pmax is the momentum 
of the most energetic electron and -Ebcam is the beam energy. 

• /Xi/yUI/ selection: Eyis/^/s < 0.7, 0.05 < PminZ-Ebeam < 0.60 and 0.05 < Pmax/-Ebeam < 

0.95, where Pmin is the momentum of the less energetic muon and Pmax is the mo- 
mentum of the most energetic muon. 

• TUTP selection: E^-^J ^ < 0.4, Pmm/-E^bcam < 0.30 and 0.04 < Pmax/-Ebeam < 0.5. 

• eufiu selection: 0.12 < E^is/V^ < 0.8 and Pmax/-E'bcam > 0.15. 

• ei/ri/ selection: Pmin/-Ebeam < 0.5 and Pmax/-Ebeam > 0.1. 

• IIUTU selection: Pram/E^eam < 0.5 and PmaxZ-^beam > 0.1. 

Table 3 shows the results for the different selections at centre-of-mass energies from 189 GeV 
to 209 GeV. Since the detection efficiencies are independent of the charged slepton species within 
the statistical uncertainties, results are quoted irrespective of the charged slepton species but 
only for the different possible final states. The lowest efficiencies arise from final states with 
two taus and missing energy and the highest ones from final states with two muons and missing 
energy. They range from 18% (TiyTiy) to 66% {^v^u) for charged slepton masses between 45 
and 103 GeV. The main contribution to the background comes from four-fermion processes 
which, depending on the final state lepton fiavours, account for 90—100% of the Standard 
Model expectation. 
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Physics process 


Eff. (%) 


Data 


Tot. bkg MC 


Data > 206 GeV 


MC > 206 GeV 


i+t 












evev 


54 - 64 


135 


145.3 ± 0.9 ± 10.7 


31 


28.6 ± 0.1 ± 3.0 




48 - 66 


99 


101.8 ± 0.6 ± 11.1 


22 


21.9 ± 0.1 ± 3.3 




18 - 27 


166 


161.5 ± 2.2 ± 22.0 


25 


35.6 ± 0.6 ± 6.7 


ev\xv 


55 - 64 


463 


502.4 ± 2.0 ± 37.1 


94 


106.2 ± 0.3 ± 10.5 


evTV 


37 - 44 


429 


459.9 ± 2.7 ± 39.7 


85 


98.0 ± 0.7 ± 11.4 




38 - 47 


434 


451.0 ± 2.7 ± 41.0 


84 


97.4 ± 0.7 d= 11.9 



Table 3: Final states with two leptons and missing energy for = 189 — 209 GeV and 
for y/s > 206 GeV: signal selection efficiencies, observed numbers of events in the data and 
background estimates from Standard Model processes. The first uncertainty on the background 
estimate is statistical and the second is systematic. 



5.3 Final states with four leptons, with or without missing energy 

Final states with four charged leptons and no missing energy may result from direct sneutrino 
decays via a A coupling while final states with missing energy may result from indirect sneutrino 
decays via a A coupling. Two analyses have been developed and optimised separately for these 
two final states. 

The following criteria were applied to select a possible signal with four leptons and missing 
energy: 

(1) The scaled event total momentum was required to be Pvis/v^ > 0.035 and the direction of 

the event missing momentum was required to satisfy | cos^missl < 0.9. 

(2) Events were required to have between three and eight tracks, each with a transverse mo- 

mentum relative to the beam axis greater than 1.5 GeV. 

(3) Events had to contain at least three identified leptons, each with a transverse momentum 

relative to the beam axis greater than 1.5 GeV. 

(4) It was also required that Ey^s/ y/s < 0.9. 

(5) The total leptonic energy, defined as the sum of the energies of all identified leptons, was 

required to be greater than 0.65 x Eyig. 

(6) It was required that the scaled longitudinal component of the event momentum satisfy 

PYis,z/Vs < 0.25. 

Table 4 shows the results for the selections with the highest and the lowest efficiencies at 
centre-of-mass energies from 189 GeV to 209 GeV. Efficiencies are quoted for Am < mp/2. 
Since the detection efficiencies are independent of the sneutrino species within the statistical 
uncertainties, results arc quoted irrespective of the sneutrino species. These analyses do not 
differentiate the difi^crcnt Icpton flavours present in the flnal states and the number of observed 
events in the data and the background estimates are therefore the same for the various signal 
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Physics process 


Eff. (%) 


Data 


Tot. bkg MC 


Data > 206 GeV 


MC > 206 GeV 


vi) — > 
liliuuiiliuv 

TTVVTTVV 


72 - 83 
13 - 28 


16 


15.6 ± 0.3 ± 2.5 


4 


4.8 ± 0.1 ± 0.7 



Table 4: Final states with four leptons and missing energy for ^/s = 189 - 209 GeV and 
separately for ^/s > 206 GeV: signal selection efficiencies for I^m < myjl, observed number of 
events in the data and background estimates from Standard Model processes. The uncertainties 
on the background estimates are the quadratic sum of the statistical and systematic uncertainties. 
Many more final states have been studied; only the final states with the highest and lowest 
efficiencies are reported here. 

final states. However, the signal detection efficiencies depend strongly on the charged lepton 
flavours present in the final states. The lowest efficiencies arise from final states with four taus 
and missing energy and the highest ones from final states with four muons and missing energy. 
They range from 13% {ttvvttvv) to 83% {ixjii/'u iijivv) for a sneutrino mass between 45 and 
103 GeV and Am < mp/2. As the value of Am is increased above mj>/2, the selection efficiencies 
drop by a factor 2 to 3, depending on m^ and Am, due to the larger fraction of missing energy 
and lower energy of the charged leptons in the final state. The main contribution to the 
background comes from four-fermion processes, which account for about 75% of the Standard 
Model expectation. The remaining 25% arises from processes yielding two-lepton final states. 

To select final states without missing energy, the following requirements were imposed: 

(1) Events were required to have at least three tracks with a transverse momentum relative to 

the beam axis greater than 1.0 GeV. 

(2) Events had to have between three and ten identified leptons, each with a transverse mo- 

mentum relative to the beam axis greater than 1.5 GeV. 

(3) It was required that -Evis/ \fs > 0.9. 

(4) To eliminate events containing photons, events were required not to have electromagnetic 

energy clusters with energies greater than 10 GeV not associated to a track. 

(5) The total leptonic energy was required to be greater than 0.7 x £^vis- 

(6) To reduce the residual four-fermion background, pairs were formed with the four most 

energetic particle tracks, and the invariant mass was computed for each pair. Events 
were selected if one of the three possible pairings satisfied \mi^j —mk^i\/ {mi^j + mk^i) < 0.4, 
where m^j is the invariant mass of the pair Only pairs with invariant mass m^j 

greater than 20 GeV were used in the computation. 

To maximise the detection efficiencies for final states without missing energy, especially 
for the four tau final states which contain missing energy from the neutrinos present in the 
tau decays, this analysis was combined with the analysis for final states with missing energy, 
previously described. Events passing either set of criteria were accepted as candidates for the 
pair-production of sneutrinos followed by a direct decay. Table 5 shows the results for the 
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Physics process 


m 


■ (%) 


Data 


Tot. bkg MC 


Data > 206 GeV 


MC > 206 GeV 




75 


- 81 
















42 


32.0 ± 0.3 ± 3.2 


8 


6.7 ± 0.1 ± 0.9 


TTTT 


32 


- 36 











Table 5: Final states with four leptons and no missing energy (including candidates from the 
analysis for four leptons with missing energy) for ^/s — 189 - 209 GeV and for > 206 GeV: 
signal selection efficiencies, observed number of events in the data and background estimates 
from Standard Model processes. The uncertainties on the background estimates are the quadratic 
sum of the statistical and systematic uncertainties. Many more final states have been studied; 
only the final states with the highest and lowest efficiencies are given in the table. 

selections with the highest and the lowest efficiencies at centre-of-mass energies from 189 GeV 
to 209 GeV. The lowest efficiencies arise for final states with four taus and the highest ones from 
final states with four muons. They range from 32% (tttt) to 81% {ii/iiiii) for a sneutrino mass 
between 45 and 103 GeV. The main contribution to the background comes from four-fermion 
processes, which account for about 80% of the Standard Model expectation. The remaining 
20% arise from processes yielding two-lepton final states. 

5.4 Final states with six leptons and missing energy 

Events with six charged leptons and missing energy in the final state may result for example 
from indirect charged slepton decays with a A coupling. Here, the strategy adopted was to 
design two complementary searches for events with six leptons either with or without miss- 
ing energy. Events satisfying either selection were accepted, which improved the detection 
efficiencies compared to only allowing events with significant missing energy. 

The following criteria were applied for the selection of final states with six leptons and 
missing energy: 

(1) The direction of the event missing momentum was required to satisfy | cos^miss] < 0.9. 

(2) The scaled event total momentum was required to be PvIs/aA > 0.035. 

(3) To reduce the background from two-photon processes and di-lepton final states, it was 

required that 0.2 < E^is/\/s < 1.0. 

(4) Events were required to have at least three tracks with a transverse momentum relative 

to the beam axis greater than 0.3 GeV. Tracks originating from photon conversions were 
excluded from this analysis. 

(5) Events had to contain between four and ten identified leptons, at least two of them with a 

transverse momentum relative to the beam axis greater than 1.5 GeV, and the third one 
with a transverse momentum relative to the beam axis greater than 0.3 GeV. 

(6) The total leptonic energy was required to be greater than 0.55 x E^^- 
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These criteria were used for the search for final states with six leptons without missing 
energy: 

(1) To reduce the background from two- photon processes and di-lepton final states, it was 

required that 0.2 < E^s.^/ ^/s < 1.2. 

(2) Events were required to have at least five tracks with a transverse momentum with respect 

to the beam axis greater than 1.0 GeV. Tracks from photon conversions were excluded. 

(3) Events had to contain between four and twelve identified leptons, each with a transverse 

momentum with respect to the beam axis greater than 1.5 GeV. 

(4) The total leptonic energy was required to be greater than 0.40 x E^is. 



Pln'sics i)r()cess 
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38 


- 63 










eimveinw 


78 
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59 


- 69 


9 


8.5 ± 0.2 ± 1.3 


1 


1.5 ± 0.04 ± 0.3 




90 


- 96 










TTTVTTTV 


24 


- 41 












79 


- 89 











Table 6: Final states with six leptons and missing energy for y/s — 189 - 209 GeV and for y/s > 
206 Ge V: signal selection efficiencies, observed number of events in the data and background 
estimates from Standard Model processes. The uncertainties on the background estimates are 
the quadratic sum of the statistical and systematic uncertainties. 



Table 6 shows the results for the events passing cither selection at centrc-of-mass energies 
from 189 GeV to 209 GeV. Since the detection efficiencies arc independent of the charged slepton 
species within the statistical uncertainties, results are quoted irrespective of the charged slepton 
species. These analyses do not differentiate the different lepton flavours present in the flnal 
states and the number of observed events in the data and the background estimates are therefore 
the same for the various signal final states. However, the signal detection efficiencies depend 
strongly on the charged lepton flavours present in the flnal states. The lowest efliciencies arise 
from flnal states with six taus and missing energy and the highest ones from flnal states with 
six muons and missing energy. They range from 24% {ttti'ttti') to 96% {miiiuiiiiiiu) for a 
charged slepton mass between 45 and 103 GeV. The background is entirely due to four-fermion 
processes. 

In the small Am region, the final state charged leptons resulting from the charged slepton 
decay into a charged lepton and a neutralino may not be detected due to the small phase space 
available. In this case, the analysis searching for four leptons and missing energy is applied 
to maximise the detection efficiencies. The lowest efficiencies in the Am < 5 GeV range were 
used. Signal selection efficiencies are of the order of 12-14%, depending on the charged slepton 
flavour. A total of 16 events was observed in the data and a background of 15.6 ± 0.3 ± 2.5 was 
estimated from Standard Model processes for y/s = 189 - 209 GeV. For ^/s > 206 GeV, the 
numbers are respectively 4 (data) and 4.8 ±0.1 ±0.7 (background). The uncertainty quoted on 
the background estimate is the quadratic sum of the statistical and systematic uncertainties. 
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5.5 Inefficiencies and systematic uncertainties for leptonic final states 

The inefficiency due to false forward detector vetoes, described in Section 5.1, caused by beam- 
related backgrounds or detector noise was estimated to range from 1.4 to 2.8%, depending on 
the data samples, from a study of randomly triggered beam crossings. The quoted efficiencies 
and background expectations take this effect into account. 

The following systematic uncertainties on the signal detection efficiencies were considered: 

1. The statistical uncertainty on the determination of the efficiency from the Monte Carlo 
simulation. 

2. The systematic uncertainty on the integrated luminosity (0.3-0.4%). 

3. The uncertainty due to the interpolation of the efficiencies (4.0%). 

4. The lepton identification uncertainty (2.4% for the muons, 3.9% for the electrons and 
4.7% for the taus). 

5. The systematic uncertainty arising from the modelling of the variables used in the selec- 
tions. This was estimated by recalculating the detection efficiencies using a new cut value 
shifted from its original value by a ratio of the means of the cut variable distributions 
of the data and of the Standard Model samples. The difference between the original 
efficiency and the new efficiency, usually less than 5%, is taken as systematic uncertainty 
due to the modelling of a cut variable. 

The systematic uncertainty due to the trigger efficiency is negligible because of the high 
lepton transverse momentum requirement. 

All cut variables are treated as independent; hence, systematic uncertainties originating 
from each variable are added in quadrature. The total systematic uncertainty was calculated 
by summing in quadrature the individual uncertainties. 



Channel 


^stat 


^ (%) 


<L (%) 


4 £ + ^„„3s OR 4 £ 
61 + E^,, OR 6 £ 


2.2 - 6.7 
1.5-8.3 
1.5 - 4.6 
0.7-5.6 


4.3 - 31.2 
7.6 - 52.7 
7.7-20.1 
6.5 - 27.8 


7.4- 18.9 
13.7 - 17.6 

7.0 - 12.7 
12.2 - 31.8 



Table 7: Relative statistical and systematic uncertainties on the selection efficiency, crf^* and 
(tI^^, and relative total uncertainty on the number of expected Standard Model background events, 
•^Nbkg' f^''^ various multi-lepton final state searches. 

The systematic uncertainty on the number of expected Standard Model background events 
is estimated in a similar way. The statistical uncertainty on this number is small, typically 
less than 1%, due to the large size of the Standard Model event samples. For final states 
with more than two leptons, an additional systematic uncertainty on the number of expected 
background events, arising from the imperfect simulation of four-fermion processes is taken into 
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account. It is determined by comparing predictions of grc4f version 2.2 including all interfering 
four-fcrmion diagrams to the default Monte Carlo predictions, which use Vermaseren for the 
two-photon processes and grc4f version 2.1 for the other four-fermion processes, and neglect the 
interference between the two. This systematic uncertainty ranges from 10 to 17%, depending 
on the analysis. 

The relative statistical and systematic uncertainties on the selection efficiency, af'^^ and 
(7^^^, and the relative total uncertainty on the number of expected Standard Model background 
events, cr^ly.^, are summarised in Table 7 for the multi-lepton final state searches presented in 
the previous sections. 

6 Final states with two jets and two leptons 

This section describes the selection for final states with two charged leptons, two jets and no 
missing energy. These final states may result from the direct decay of pair-produced stops via a 
A coupling. In contrast to the purely leptonic final states described in the previous section, the 
topologies searched for in this analysis involve hadronic jets; more stringent cuts are needed 
to obtain a purer lepton sample. Particles are considered as electrons or muons if they are 
either identified by the selection algorithms described in [49] and [50] , or by an algorithm used 
for selecting semileptonic W decays, as described in [52]. A Neural Net (NN) based on track 
properties [53] is used to identify taus. Events were preselected by requiring the following 
criteria to be satisfied (the same criteria were also used for the analysis presented in Section 7): 

The ratio of good tracks satisfying the criteria of Section 3 to the total number of recon- 
structed tracks had to be greater than 0.2 to reduce beam-gas and beam-wall background 
events. Events with fewer than eight good tracks were not considered, in order to reduce the 
background from Bhabha scattering. Events had to contain at least one identified electron or 
muon with a momentum p greater than 3 GeV, to reduce the background from final states with 
low energy leptons (electron or muon) arising for instance from semi-leptonic quark decays or 
misidentified particles. To reduce background from two-photon processes, it was required that 
the visible energy normahsed to the centre-of-mass energy £'vis/-\/s be greater than 0.3. 

The following cuts were then apphed: 

(1) The visible energy had to be close to the centre-of-mass energy, i.e. 0.75 < £'vis/-\/s < 1.25 

for electron and muon and 0.5 < £'vis/-\/i < 1.0 for tau final states. 

(2) It was required that four jets be reconstructed using the Durham [54] algorithm, with 

1/34 > 0.001, where 1/34 is the jet resolution at which the number of jets changes from 3 to 
4. Both hadronic and leptonic objects are used in the jet reconstruction. 

(3) Events had to contain at least one pair of identified oppositely-charged lepton candidates 

of the same flavour. 

(4) To make use of the signal topology of two leptons and two jets, where a lepton and a jet 

stem from the same object, a five-constraint (5C) kinematic fit was performed for the two 
possible combinations of each lepton with each jet. The kinematic constraints are: the 
vector sum of all momenta has to be equal to zero, the total energy of all objects has to 
be equal to the centre-of-mass energy and the masses of the two reconstructed particles 
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have to be equal. From the three most energetic leptons of the same flavour, the two 

most isolated^ were selected and the rest of the event was reconstructed as two jets. The 
combination with the highest fit probability was selected. The probability for the fit, 
based on the x^, was required to be larger than 0.01. 

(5) The scaled momenta p/y^ of the most and second most energetic leptons had to be greater 

than 0.082 and 0.055, respectively, for final states with electrons and muons and greater 
than 0.055 and 0.0275, respectively, for taus. 

(6) It was required that there be no track within a cone of half opening angle of 15° around 

the most energetic lepton candidate. 

Figure 2 shows the distribution of the reconstructed stop mass after cut (6). The kinematic 
fit achieves a very good mass resolution. For stop masses above 75 GeV, the mass resolution 
is better than 0.8 GeV for final states containing either electrons or muons and better than 
1.3 GeV for final states with taus. Events arc counted only in a two sigma mass window around 
the mean reconstructed stop mass. For masses near the kinematic limit the mass resolution 
approaches zero. If the mass resolution is below 0.5 GeV, it is set to 0.5 GeV. 




Reconstructed mass (GeV) 

Figure 2: Final states with two jets and two leptons: distribution of the reconstructed stop mass 
after cut (6) for the direct decay of stops into final states with two jets and two electrons. The 
signal is shown as the dashed line for a stop mass of 75 Ge V. Data are shown as points and 

the sum of all Monte Carlo background processes is shown as a solid line. Data used for this 
plot were collected at centre- of-mass energies of 189-209 GeV. The simulated signal events have 
arbitrary normalisation. Uncertainties are statistical only. 

^The most isolated lepton is the one with the largest angle to the closest other track. 
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These cuts, including the cut of the two sigma mass window, yield an efficiency of more than 
46% for a stop mass above 75 GeV for final states with either electrons or muons. For final states 
with taus, efficiencies above 15 % are reached for stop masses above 75 GeV. Table 8 shows the 
efficiencies, the number of candidate events observed in the data and the number of background 
events for the selection at centre-of-mass energies from 189 GeV to 209 GeV, corresponding 
to an integrated luminosity of 599.6 pb~^. The main contribution to the background comes 
from four-fermion processes which amount to 100% of the Standard Model expectation for final 
states with muons and to 75—95% for final states with taus. In the electron channel, four- 
fermion events account for 45—75% of the Standard Model background, two-fermion processes 
contribute less than 5% and the remaining events are expected from two-photon interactions. 



Physics process 


Eff. (%) 


Data 


Tot. bkg MC 


Data > 206 GeV 


MC > 206 GeV 


it 

2 jets + 6^+6" 

2 jets + /i^ + /i^ 
2 jets + T+ +T- 


46 - 55 

47 - 55 
15 - 25 


13 
7 

30 


13.8 ± 1.4 ± 1.2 
7.1 ±0.4 ± 1.5 
34.4 ± 1.0 ±4.1 


3 

7 


3.3 ±0.5 ±0.3 
1.5 ±0.2 ±0.1 
7.9 ±0.6 ±0.9 



Table 8: Final states with two jets and two leptons for i/s = 189 - 209 GeV and for y/s > 
206 GeV: signal selection efficiencies, observed number of events in the data and background 
estimates from Standard Model processes. The first uncertainty on the background estimate is 
statistical and the second is systematic. The efficiencies are given for stop masses greater than 
75 GeV for the signal at ^ = 207 GeV. 

Systematic uncertainties 

The following systematic uncertainties on the signal efficiencies were considered: 

1. The statistical uncertainties from the hmited size of the Monte Carlo samples. 

2. The systematic uncertainties on the integrated luminosity (0.3 - 0.4%). 

3. The uncertainty due to the interpolation of efficiencies for mass values between the gen- 
erated stop masses (less than 3%). 

4. For the lepton identification, the highest uncertainty (4% for electrons) was taken for all 
fiavours. 

5. The systematic uncertainty arising from the modelling of the variables used in the selec- 
tion, as described below. 

6. The fragmentation of the stop was simulated using the fragmentation function from Pe- 
terson et al. with the e parameter extrapolated from measurements of charm and bot- 
tom [55]. To check the model dependence of the fragmentation, it was also performed 
using the function from Bowler [56]. No significant change in the efficiency due to the 
difference in the fragmentation function was found. The difference was at most 0.5%, 
where a variation of the e parameter of the t in the Peterson et al. scheme was included. 
This uncertainty on was propagated from the uncertainty of and the uncertainty on 
the b-quark mass as described in detail in Ref. [55]. 
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7. The signal events were produced for a zero mixing angle between the two stop eigenstates. 
The mixing angle describes the coupling between the stop and the Z°, and therefore the 
energy distribution of the initial state radiation depends on this mixing angle. To check 
the dependence of the detection efficiency on this angle, events were generated with 
9i =0.98 rad, where the stop decouples from the Z°. The change in efficiency was less 
than 0.5% for the two extreme cases. 

8. The Fermi motion of the spectator quark in the stop-hadron influences its measured mass. 
The Fermi motion was increased from 220 MeV to 520 MeV and the efficiency changes 
by no more than 1%, which was taken as a systematic uncertainty. 

The systematic uncertainty due to the uncertainty in the trigger efficiency was estimated 
to be negligible, because of the requirement of at least eight good tracks. 

All of the above contributions to the systematic uncertainties result in a total systematic 
uncertainty of less than 7%. 

The uncertainty in the signal efficiency and in the number of background events due to the 
modelling of the selection variables was estimated in the following way. First the mean and 
the r.m.s. of the distribution of the variable (y) was determined for the data and the Standard 
Model background simulation (SM MC). Then for each MC event (both signal and background) 
the variable was transformed as 

y' ^{y- mean(SM MC)) • r.m.s. (data) /r.m.s. (SM MC) + mean(data) 

The standard analysis was apphed to signal and Standard Model background events with all 
variables transformed; the difference in the MC expectation due to the transformation was 
considered as the systematic uncertainty. The contributions for the different variables were 
added in quadrature. 

The systematic uncertainty on the expected number of background events estimated in 
this way was found to be less than 20%, the largest contribution (~ 15%) stemming from the 
variable ^34. In addition, the difference in the number of selected events by comparing different 
Monte Carlo generators, as described in Section 4, was also taken as a systematic uncertainty. 

7 Final states with four jets and at least two charged 
leptons 

This section describes the event selection for final states resulting from the indirect decay of 
selectrons, smuons and staus via the coupling A . The final states consist of two charged leptons 
of the same flavour as the sleptons plus the decay products of the two neutrahnos. These will 
be two jets plus a neutral or charged lepton for each neutralino. The identification of leptons 
and the preselection are the same as described in Section 6. The selection cuts were as follows: 

(1) Lower and upper cuts in the ranges 0.5-0.7 and 0.75-1.2, respectively, were applied to the 
visible energy scaled by the centre-of-mass energy Eyis/y/s, depending on the expected 
number of neutrinos. In addition, if some missing momentum was expected, a cut was 
made on the angle of the missing momentum with respect to the beam direction at 
I cos^l < 0.95. 
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(2) The jets in the event were reconstructed using the Durham algorithm. The jet resolution 

parameter ^45 at which the number of jets changes from four to five was required to be 
greater than 0.002. The log(y45) distribution is shown in Figure 3. It demonstrates the 
good discriminating power of this quantity. The visible disagreement in the peak position 
between data and Standard Model expectation was taken into account as a systematic 
uncertainty. This cut takes into account the high multiplicity of the signal events. 

(3) At least two leptons of the flavour of the slepton had to be identified. To retain sensitivity to 

small mass differences between the slepton and the neutralino, the required scaled energy 
greater than 0.022 for the two electrons in the selectron case, the required momentum had 
to be greater than 4 GeV for both muons in the smuon case and the required momentum 
had to be greater than 3 GeV for both taus in the tau case, respectively. 

(4) In addition to the leptons required in the previous cut (3), the leptons from the neutralino 

decay had to be identified. If two additional charged leptons with a different flavour to 
the slepton were expected, both had to be identifled except in the case of two taus, where 
only one had to be identifled. If a total of four leptons of the same flavour was expected, 
including those in cut (3), only three of them had to be identified. If only one additional 
lepton was expected, it had to be identified. 

The scaled energy or momentum of the most energetic lepton had to be above a cut value 
varying between 0.044 and 0.08, depending on the topology. If a total of four leptons was 
required, for the second most energetic lepton a scaled energy or momentum larger than 
a cut value varying between 0.016 and 0.022, depending on the topology, was required. 

(5) To make use of the isolation of the leptons in the signal, one or two of the identified leptons, 

depending on the expected topology, must be isolated. The isolation criterion was that 
there be no charged track within a cone of half opening angle 93, such that | cos 99 1 = 0.99, 
around the track of the lepton. 

Table 9 shows the efficiencies, the number of candidate events observed in the data and 
the number of background events for the selection at centre-of-mass energies from 189 GeV to 
209 GeV, corresponding to an integrated luminosity of 599.6 pb~^. Good agreement was found 
between the numbers of events expected and observed. The largest difference is observed for 
the final state with two taus and two neutrinos, where there is a 1.8 cXstat excess in the number 
of data events. The main contribution to the Standard Model background comes from four- 
fermion processes; two-fermion multi-hadron events contribute up to 30% and other processes 
were negligible. 

Systematic uncertainties 

The systematic uncertainties on the signal efficiencies and on the expected number of back- 
ground events were estimated in the same manner as in Section 6: 

1. The statistical uncertainty from the hmited size of the Monte Carlo samples. 

2. The systematic uncertainty on the integrated luminosity (0.3 - 0.4%). 



23 




io8(y45) 

Figure 3: Final states with four jets and at least two charged leptons: distribution of the variable 
^45 before the cut on this variable for the indirect decay of selectrons into final states with four 
jets and four electrons. The signal is shown as the dashed line for a selectron mass of 102 GeV 
and a mass difference Am — m/2. Data are shown as points and the sum of all Monte Carlo 
background processes is shown as a solid line. Data used for this plot were collected in the 
year 2000 at a luminosity weighted average centre- of-mass energy of 206.6 GeV. The visible 
disagreement in the peak position between data and Standard Model expectation is taken into 
account as a systematic uncertainty. The simulated signal events have arbitrary normalisation. 
The arrow points into the region accepted by the cut. Uncertainties are statistical only. 
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Physics process 


Eff. (%) 


Data 


Tot. bkg MC 


Data > 206 GeV 


MC > 206 GeV 


e^e" — > 
4 jets + 2e + 2e 
4 jets + 2e + 2;U 
4 jets + 2e + 2r 
4 jets + 2e + ez/ 
4 jets + 2e + /iz/ 
4 jets + 2e + Ti/ 
4 jets + 2r + 2?/ 


39 - 69 

25 - 53 
6-25 
17 - 40 
21 - 52 
2-12 
6 - 17 


13 
2 

14 
16 
16 
7 

6 


13.3 ± 0.6 ± 1.2 

2.6 ±0.3 ±0.4 
12.7 ± 0.6 ± 1.6 
17.6 ± 0.7 ± 1.5 
12.3 ± 0.6 ± 1.4 
7.3 ±0.5 ±0.9 
5.8 ±0.1 ±0.1 


2 

5 
2 
4 
1 
2 


2.5 ±0.3 ±0.3 
0.7±0.1 ±0.1 

2.7 ±0.3 ±0.2 
3.3 ±0.3 ±0.3 

2.8 ±0.3 ±0.4 
1.5 ±0.2 ±0.2 
1.1 ±0.3 ±0.1 


— > 

4 jets + 2^ + 2e 
4 jets + 2/i + 2/i 
4 jets + 2/x + 2t 
4 jets + 2/i + ei/ 
4 jets + 2ii + iiv 
4 jets + 2/x + rz/ 
4jets + 2^ + 2z/ 


34 - 54 
52 - 79 
9-32 
25 - 52 
30 - 51 
6-22 
15 - 35 


2 
4 
2 
5 
4 
2 
6 


1.6 ±0.2 ±0.2 
2.9 ±0.3 ±0.6 
3.6 ±0.3 ±0.5 
4.4 ±0.3 ±0.8 
3.6 ±0.3 ±0.6 
2.3 ±0.2 ±0.3 
3.8 ±0.3 ±0.4 



2 


2 




0.2 ±0.1 ±0.0 
0.7 ±0.1 ±0.2 
0.7±0.1±0.1 
1.0 ±0.2 ±0.3 
0.9 ±0.2 ±0.2 
0.4 ±0.1 ±0.0 
0.8 ±0.1 ±0.0 


f^f^ — > 
4 jets + 2r + 2e 
4 jets + 2r + 2^ 

4 jets + 2r + 2r 
4 jets + 2r + ez/ 
4 jets + 2r + /xz/ 
4 jets + 2r + rz/ 

4 jets + 2- + 2// 


19 - 52 
19 - 56 

7- 22 
4-16 
4-14 

8- 21 

3 - 8 


9 
7 
53 
15 
12 
26 
83 


9.4 ±0.5 ±0.6 
8.8 ±0.4 ±1.4 

46.0 ± 1.1 ±3.9 

12.1 ±0.6 ±1.1 

12.8 ±0.5 ±1.4 

24.9 ±0.8 ±2.3 

66.3 ± 1.2 ± 1.9 


1 


6 
1 
1 

6 

16 


1.9 ±0.3 ±0.1 

2.2 ±0.2 ±0.5 

9.3 ±0.6 ±0.6 
2.1 ±0.3 ±0.2 
2.8 ±0.3 ±0.4 
5.1 ±0.4 ±0.5 
13.6 ± 0.6 ± 0.2 



Table 9: Final states with four jets and at least two charged leptons for y/s — 189 - 209 GeV 
and for y/s > 206 GeV. Signal selection efficiencies, observed number of events in the data and 

background estimates from Standard Model processes. The first uncertainty on the background 
estimate is statistical and the second is systematic. The efficiencies are given for slepton masses 
larger than 45 GeV and for a signal at y/s = 207 GeV. 
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3. The uncertainty due to the interpolation of efficiencies for mass values between the gen- 
erated mass points (4%). 

4. For the lepton identification, the highest uncertainty (4% for electrons) was taken for all 
flavours. 

5. The systematic uncertainty arising from the modelling of the variables used in the selec- 
tion. This results in a total systematic uncertainty on the efficiency of 5.2 - 13.4% for the 
selectron selection, of 5.2 - 12.0% for the smuon selection and of 5.2 - 19.3% for the stau 
selection. 

6. Prom the studies of the fragmentation in Section 6, the systematic uncertainty for this 
analysis was estimated to be less than 1%. 

The systematic uncertainty due to the uncertainty in the trigger efficiency was negligible, 
because of the requirement of at least eight good tracks. 

The systematic uncertainty on the expected number of background events was again esti- 
mated to be less than 20% for all cases. 



8 Final states with four jets and missing energy 

Indirect decays of sneutrinos via a A' coupling can lead to flnal states with four jets and a large 
amount of missing energy due to the four undetected neutrinos. The dominant background 
contribution comes from four-fermion processes, mainly W+W^ qqiu, and radiative or mis- 
measTired two-fermion events. The selection procedure closely follows the one described in [16] 
with small changes required at the higher centre-of-mass energies. The selection procedure is 
described below: 

(1) The event had to be classified as a multi-hadron final-state as described in [57]. 

(2) The visible energy of the event was required to be less than 0.75 x ^/s to account for the 

undetectable neutrinos in the final state. 

(3) To reject two-photon and radiative two-fermion events, the transverse momentum of the 

event should be larger than 0.075 x ^/s, the total energy measured in the forward calorime- 
ter, gamma-catcher and silicon tungsten calorimeter should be less than 15 GeV, and there 
should be no significant energy deposit in the scintillating tile counters. The missing mo- 
mentum should not point along the beam direction ([ cos^missl < 0.96). 

(4) The events were forced into four jets using the Durham jet-finding algorithm, and rejected 

if the jet resolution parameter 1/34 was less than 0.0008. All jets must contain at least one 
charged particle track. 

(5) An additional cut was applied against semi-leptonic four-fermion events, vetoing events 

containing isolated leptons. The lepton veto is based on a NN routine [53], which was 
designed to identify tau leptons. The algorithm to select 1-prong tau candidates looks 
for high momentum isolated tracks, and therefore it is also suitable for vetoing leptons of 
other flavours. If any lepton candidate was found with a NN output larger than 0.97, the 
event was rejected. 
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Finally, a likelihood selection was employed to classify the remaining events as two-fermion, 
four-fcrmion or signal-like processes. The method is described in [58]. The signal reference 
histograms were produced separately for first generation sneutrinos where t-channel produc- 
tion is also expected and for second or third generation sneutrinos with s-channel production 
only. For a given centre-of-mass energy, all the generated signal events with the coupling A^2i > 
whatever their masses, were used with equal weight to form the reference distributions. The 
information from the following variables was combined: 

• the effective centre-of-mass energy [59] of the event; 

• the transverse momentum of the event; 

• the cosine of the polar angle of the missing momentum vector; 

• the D event-shape parameter [60] which provides information about the planar nature of 
an event; 

• the aplanarity of the event which measures the transverse momentum component out of 
the event plane and which is defined as 3/2 times the smallest eigenvalue of the sphericity 
tensor; 

• the logarithm of ^34; 

• the highest track momentum; 

• the highest electromagnetic cluster energy; 

• the number of lepton candidates in the event with a NN output larger than 0.5; 

• the mass of the hadronically decaying W after a kinematic fit to the W''"W~ — > qq£u 
hypothesis; 

• the cosine of the smallest jet opening angle, defined by the half-angle of the smallest cone 

containing 68% of the jet energy; 
The distribution of the likelihood output is shown in Figure 4 for the data, the estimated 
background and a simulated signal sample in the electron sneutrino selection. The event was 
rejected if its likehhood output was smaller than 0.95. 

Table 10 shows the results of the different selections on the data collected at centre-of-mass 
energies between 189 GeV and 209 GeV corresponding to an integrated luminosity of 610.6 pb~^. 

The efficiencies shown in Table 10 arc quoted for the coupling Am < mo/2. For 

Am > mp/2, the detection efficiency drops by a factor 2 to 4, depending on ma and Am, due 
to the larger fraction of missing energy and lower energy final state jets. For a mass difference of 
5 GeV, the primary neutrinos carry away less energy and the efficiency increases by an absolute 
value of about 10%. The efficiency is small for fight sneutrino masses due to the effect of the 
initial-state radiation and the larger boost of the particles, which make the event similar to the 
QCD two-fermion background. About 94% of the events expected from the Standard Model 
originate from four-fermion processes. 



Physics process 


Eff. (%) 


Data 


Tot. bkg MC 


Data > 206 GeV 


MC > 206 GeV 


i>ei>e 4 jets + Emiss 


1 - 25 


54 


39.5 ± 0.5 ± 4.3 


8 


8.2 ± 0.2 ± 1.3 


VtVt 4 jets + Emiss 


1 - 22 


57 


39.4 ± 0.5 ± 4.3 


9 


8.3 ± 0.2 ± 2.0 



Table 10: Final states with four jets and missing energy for ^/s = 189 - 209 GeV and for ^/s > 
206 GeV: signal selection efficiencies for sparticle masses above 45 GeV and Am < mo/2, 
observed numbers of events in the data and background estimates from Standard Model processes. 
The first uncertainty on the background estimate is statistical and the second is systematic. The 
numbers of observed events are highly correlated between the selections. 
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Figure 4: Final states with four jets and missing energy, electron sneutrino selection at ^/s = 
189-209 GeV: distribution of the likelihood output for the selected events in the data (points with 
error bars), the estimated background normalised to the integrated luminosity of the data (solid 
line) and a simulated sample of electron sneutrino pair-production for m^^ — 102.5 GeV with 
arbitrary normalisation (dashed line). Uncertainties are statistical only. The arrow indicates 
the cut position. 
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The observed approximately 2 Ugtat excess in the data spreads over the different centre-of- 
mass energies and it is concentrated in the 0.075 < px/v^ < 0.15 region, not favoured by 
events with pair-produced high mass sleptons. An excess appears from the first steps of the 
selection. After cut (5), 1445 events were observed for 1240 expected from Standard Model 
processes. As seen in Figure 4, there is a significant excess in the low likelihood region, which 
is expected to be almost equally populated by two- and four-fermion events. 

Systematic uncertainties 

The following sources of uncertainties were considered on the selection efficiencies (all quoted 
uncertainties are relative): 

1. The limited MC statistics (5.4-28.7%). 

2. The statistical and systematic uncertainties on the luminosity measurement (0.3—0.4%). 

3. The modeUing of the preselection variables (2.2—2.5%). 

4. The modelling of the variables used in the likelihood selection (20.9—23.4%). 

5. The modelling of the lepton veto (1%). 

6. The interpolation of efficiencies (2.7—5.2%). 

Similarly the background estimate is affected by uncertainties due to: 

1. The hmited MC statistics (1.2%). 

2. The statistical and systematic uncertainties on the luminosity measurement (0.3—0.4%). 

3. The modelling of the preselection variables (8.6—9.6%). 

4. The modeUing of the variables used in the hkelihood selection (27.8—29.1%). 

5. The modeUing of the lepton veto (1%). 

6. The imperfect knowledge of Standard Model cross-sections (2%). 

7. The imperfect simulation of four-fermion processes determined by comparing the predic- 
tions of grc4f, KORALW and KandY (16%). 

The total systematic uncertainty excluding MC statistics ranges from 21.2 to 24.3% on the 
signal detection efficiency and from 33.3 to 34.6% on the background rate. 

The effect of the modelling of the preselection and likelihood variables was estimated for 
each variable using the transformation described in Section 6. Applying the same analysis 
with unchanged reference histograms to calculate the hkehhood, the difference in the MC 
expectation due to the transformation was considered as the systematic uncertainty. The 
largest contributions both for the signal and for the background come from the modelling of 
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the logarithm of ^34, the energy of the most energetic electromagnetic cluster and the cosine of 
the smallest jet opening angle in the likehhood selection. 

The uncertainty introduced by interpolating the signal efficiencies between simulated signal 
masses was estimated in the following way: the signal efficiencies were randomly smeared 
within their statistical uncertainty and the resulting numbers were used in the interpolation. 
This process was repeated 100 times, and the absolute differences between the original and the 
smeared fits were averaged. For each slepton mass, the mean difference was taken as a measure 
of the systematic uncertainty. There is no interpolation between the generated centre-of-mass 
energies: the efficiency determined from the signal sample with closest generated centre-of- 
mass energy was always taken. The largest difference between the actual and the generated 
centre-of-mass energies is about 3 GcV, and its effect on the result is negligible. 

The inefficiency due to the forward energy veto, by which both the signal efficiencies and 
the background estimates were decreased, was found to vary between 1.9 and 3.7% by studying 
randomly triggered events. 

9 Final states with four jets and no missing energy 

Direct decays of sleptons via A' couplings and squarks via A" couplings can result in final states 
with four well-separated, high multiplicity hadron jets and a large amount of visible energy. 
The main background comes from four-fermion processes, dominantly W'*"W~ — > qqqq, with 
some contribution from qq(7) events with hard gluon emission. 

The analysis is similar to the one of [16]. It consists of a cut-based preselection to reduce the 
two-photon and two-fermion backgrounds and a likelihood selection to suppress the contribution 
from Standard Model four-fermion events. The preselection consists of the following steps: 

(1) The event has to be classified as a multi-hadron final state as described in [57]. 

(2) The effective centre-of-mass energy of the event, [59], was required to be greater than 

0.82 X to reject events with large initial state radiation. 

(3) To ensure that the events were well contained in the active region of the detector, the 

visible energy should be larger than 0.7 x ^/s. 

(4) The events were forced into four jets using the Durham jct-fiiiding algorithm, and rejected 

if the jet resolution parameter ^34 was less than 0.0025. Moreover, all jets must contain 
at least one charged particle. 

(5) A four-constraint kinematic fit (4C-fit) requiring energy and momentum conservation for 

the jet four- momenta should yield a probability larger than 10~^. 

(6) To test the compatibility with pair-produced equal mass objects and to obtain the best 

possible resolution for the reconstructed mass of a hypothetical sfermion, the jet four- 
momenta were refitted requiring energy and momentum conservation and equal jet pair 
invariant masses (5C-fit). The event was kept if at least one of the three jet pairing 
combinations has a probability larger than 10~^. The combination with the highest 
probability was considered later in the mass reconstruction. 
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(7) A cut was applied on the C event shape parameter [60] which provides an effective measure 
of the multi-jet structure of the event, C > 0.45. 

Finally, a likelihood selection was employed to classify the remaining events. Three event 
classes were used: signal, two-fermion and four-fermion. The signal reference histograms were 
produced separately for selectrons where i-channel production plays an important role, for 
electron-sneutrinos where the t-channel process also contributes but to a lesser extent and for 
second or third generation neutral and charged slcptons with s-channel production only. Since 
squarks are coloured particles they were treated separately. For a given centre-of-mass energy, 
all the generated masses were used with equal weight to form the reference distributions. The 
following variables were used as inputs to the likelihood calculation: 

• the cosine of the polar angle of the thrust axis; 

• the cosine of the smallest angle between the directions of any two of the four reconstructed 

jets; 

• the difference between the largest and smallest jet energy after the 4C-fit; 

• the smallest difference between the reconstructed masses of the two jet pairs from any of 
the three possible jet pair combinations; 

• the cosine of the direction of the jet pair momentum multiplied by the charge of the jet 
pair^ for the combination with the highest probability given by the 5C-fit. 

Events were accepted if their likelihood output was larger than 0.5 for selectrons, 0.55 for 
electron sneutrinos, 0.6 for second and third generation neutral and charged sleptons and 0.7 
for squarks. 

Table 11 shows the numbers of selected data and expected background events for the dif- 
ferent selections at centre-of-mass energies from 189 GeV to 209 GeV, corresponding to an 
integrated luminosity of 610.6 pb~^. More than 92% of the events expected from the Standard 
Model originate from four-fermion processes. Figure 5 shows, as an example, the jet pair in- 
variant mass distribution of the selected events after the 5C-fit in the selectron selection: the 
dominance of W+W~ production in the Standard Model background is clearly visible from the 
strongly peaked distribution around the W boson mass. The jet pair invariant mass resolution 
for signal events is o" ~ 0.6 — 1.6 GeV, depending on the type and the mass of the produced spar- 
ticle. To profit from the good mass resolution, events were selected if they were in a ±2(T mass 
window around the test mass. The efficiencies in Table 11 are given after the mass selection 
and all correspond to the Yukawa coupling giving the worst result. 

Systematic uncertainties 

The systematic uncertainties were estimated in the same manner as in the search for four jets 
and missing energy. For the selection efficiency the following sources of uncertainties were 
considered: 

1. The limited MC statistics (3.1-11.9%). 

2. The statistical and systematic uncertainties on the luminosity measurement (0.3—0.4%). 

^The charge of the jet pair was calculated as I]g(j)PL(i) > where the sum goes over each track within the two 
jets, q{i) is the charge of the track and pL(i) is its momentum parallel to the jet direction. A charge of was 
assigned to the jet pair with the larger charge, and a charge of —1 to the other. 
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Physics process 


Eff. (%) 


Data 


Tot. bkg MC 


Data > 206 GeV 


MC > 206 GeV 


e'^e~ — > 4 jets 


14 - 46 


922 


917 ± 10 ± 69 


193 


210 ± 5 ± 38 


i>e^'e 4 jets 


15 - 43 


760 


762 ± 8 ± 56 


158 


175 ± 4 ± 30 














i>ijp^, UrU^ 4 jets 


9-39 


584 


589 ± 6 ± 43 


123 


138 ± 3 ± 24 


qq ^ 4 jets 


9-37 


393 


393 ± 5 ± 29 


81 


89 ± 2 ± 15 



Table 11: Final states with four-jets and no missing energy for ^/s — 189 - 209 GeV and for 
a/s > 206 GeV: signal selection efficiencies for sparticle masses above 45 GeV, observed num- 
bers of events in the data and background estimates from Standard Model processes. The first 
uncertainty on the background estimate is statistical and the second is systematic. The efficien- 
cies are given after the mass selection described in the text while the data and the background 
correspond to the whole mass region. The inefficiency related to the fragmentation model was 
taken into account for the squark result. There is a significant correlation amongst the different 
selections at the same centre- of-mass energy. 



200 I ' I ' ' I ' I I ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' 




Jet pair invariant mass (GeV) 

Figure 5: Final states with four jets and no missing energy, selectron selection at y/s = 189- 
209 GeV: distribution of the reconstructed mass for the selected events in the data (points with 
error bars), the estimated background normalised to the integrated luminosity of the data (solid 
line) and a simulated sample of selectron pair-production for mg = 62 GeV with arbitrary 
normalisation (dashed line). Uncertainties are statistical only. The arrows indicate the mass 
window which would be selected for the specific case of the signal shown. 
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3. The modelling of the preselection variables (15.9—18.2%, largest contribution of 10.2—15.3% 
from the probability of the 4C-fit). 

4. The modeUing of the variables used in the hkelihood selection (6.6—10.9%). 

5. The interpolation of efficiencies (1.6—3.2%). 

The background estimate is subject to uncertainties due to 

1. The limited MC statistics (1.1%). 

2. The statistical and systematic uncertainties on the luminosity measurement (0.3—0.4%). 

3. The modelling of the preselection variables (13.2—16.3%, largest contribution of 8.7— 13.3% 
from the probability of the 4C-fit). 

4. The modeUing of the variables used in the hkelihood selection (4.5—10.7%). 

5. The imperfect knowledge of Standard Model cross-sections (2%). 

6. The imperfect simulation of four-fermion processes (3.3—6.4%). 

The total systematic uncertainty excluding MC statistics ranges from 17.3 to 21.5% on the 
signal detection efficiency and from 16.4 to 18.3% on the background rate. 

As explained in Section 4, the squark pair Monte Carlo samples were produced using an 
independent fragmentation model which gives narrower, better separated jets and therefore 
optimistic selection efficiencies. To estimate this effect, the efficiencies were compared for 
smuon Monte Carlo samples with independent and string fragmentation models. From this 
study an additional relative inefficiency of 7.9% to 36.9% was derived, by which the squark 
detection efficiencies were decreased. For the same reason, the squark results at 183 GeV [16] 
were updated using a similar procedure and these new results are used later in the combination. 

10 Interpretation 

No significant excess of signal-like events was observed in the data with respect to the expected 
background for any analysis listed in Table 1. Production cross-section and mass hmits were 
therefore computed using the data collected at LEP from centre-of-mass energies (i/i) between 
189 GeV and 209 GcV and then combined with results obtained using data at lower energies [16]. 
These limits take into account indirect limits obtained from the study of the Z° width at LEPl 
and therefore concern only sparticle masses above 45 GeV. All limits presented here are quoted 
at the 95% confidence level. 

Two approaches were used to present sfermion production limits. In the first, upper limits on 
the production cross-sections were calculated as functions of the sfermion masses with minimal 
model assumptions. These upper hmits in general do not depend on the details of the SUSY 
models, except for the assumptions that the sparticles are pair-produced and that only one 
A-like coupling at a time is non-zero [16]. In the second approach, limits on the sfermion 
masses were calculated in an ii-parity violating framework analogous to the CMSSM, where 
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mass limits were derived for tan/3 = 1.5 and ji — —200 GeV. This choice of parameters, which 
generally results in small sfcrmion production cross-sections, is a convenient benchmark for 
limit setting and is used by other collaborations, so results may be compared and combined. 
For the indirect sfermion decays, we used the branching ratios for the decay / — > /x5 predicted 
by the CMSSM, and we conservatively assumed no experimental sensitivity to any other decay 
mode. The branching ratio for direct decay was always assumed to be unity, as only one A 
coupling at a time was allowed to be different from zero. 

As in [16], the relative branching ratios of the neutralino into a final state with a charged 
or a neutral lepton were varied between and 1 to avoid a dependence of the results on the 
CMSSM parameters. A likelihood ratio method [61] was used to determine an upper limit for 
the cross-section. This results in a cross-section limit as a function of the branching ratio and 
the sfermion mass. A limit independent of the branching ratio was determined by taking the 
lowest limit at each sfermion mass. For the direct decays, the final states are fully determined 
by the indices of the coupling considered. 

In the following sections, cross-section limits are shown for the various direct and indirect 
decays studied in this paper and are listed in Table 1. In each cross-section plot for the decays 
via A couplings, the curves corresponding to the highest and the lowest cross-section limits are 
shown. Generally, the highest excluded cross-section comes from final states with a maximum 
number of muons and no taus, while the worst results come from final states with many taus, 
due to their lower detection efficiency. In the other cases, only the curve corresponding to the 
worst cross-section limit is shown amongst all possible cross-section limits resulting from the 
couplings considered. 

In the CMSSM framework, the exclusion regions for the indirect decays are valid for m^o > 
10 GeV, to ensure prompt decays. The region in the plane (m^, m^o) corresponding to m^o < 
10 GeV, where the hfetime of sparticles would be sufficiently long to produce a secondary 
decay vertex, clearly detached from the primary vertex, or even outside the detector, is labelled 
"Lifetime signature" and is not excluded by these searches. 

For indirect sparticle decays, mass lower limits are quoted for Xi = 10 GeV, referred to as 
a low-mass Xi- In the small Am region (Am = m^ — m^o < 5 GeV), the final state charged 
leptons resulting from the charged slepton decay into a charged lepton and a neutralino may 
not be detected due to the small phase space available. This region is excluded for the charged 
sleptons decaying indirectly via a A coupling. For the charged sleptons decaying indirectly via a 
A coupling, no exclusion was possible in this region. The exclusion region for the direct decays 
is independent of Am. 

The production cross-section for left-handed charged sleptons is always larger than that for 
right-handed charged sleptons; therefore, whenever applicable for the superpotential form given 
in Equation (1), results are conservatively quoted for right-handed charged sleptons. 

10.1 Selectron limits 

Upper limits on the cross-sections of pair-produced selectrons followed by the direct decay via 
a A coupling are shown in Figure 6. Due to the structure of the LLE term and the presence of 
a neutrino in the decay products, only left-handed sleptons can decay via this term, the other 
decays being suppressed. Since the selection efficiency for this analysis only depends on the 
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final states and not on the parent slepton type, this limit also holds for the pair production of 
smuons and staus decaying directly via a A coupling. 

OPAL 

£ 0.4 
^ 0.3 

0.2 

0.1 

^ 50 60 70 80 90 100 

m(\^) (GeV) 

Figure 6: Charged slepton direct decays via a X coupling: upper limits at the 95% C.L. on the 
pair-production cross- sections of left-handed sleptons decaying directly. The dashed line shows 
the lowest upper limit (tutu final states) while the solid line shows the highest one (fiufiu final 
states). 

Upper limits on the cross-sections of pair-produced selectrons decaying indirectly via a A 
coupling are shown in Figure 7. 

Upper limits on the cross- sect ions of pair-produced selectrons decaying directly via a A' 
coupling to a four-jet final state are shown in Figure 8. The peak structure visible in the 
figure at approximately the mass of the W-boson comes from irreducible background due to W 
pair-production. 

Upper limits on the cross-sections of pair-produced selectrons decaying indirectly via a A' 
coupling are shown in Figure 9 for the indirect decay of a e^ in the electron channel, in the 
muon channel and in the tau channel. 

In the CMSSM, the selectron pair-production cross-section is enhanced by the presence of 
a t-channel diagram. The exclusion limits for selectrons decaying via a A coupling are shown 
in Figure 10(a). In the region where the neutralino is heavier than the selectron, only direct 
decays are possible. When the neutralino is lighter than the selectron, the indirect decays are 
expected to be dominant. The exclusion refers to right-handed selectrons for the indirect decays 
and to left-handed selectrons for direct decays. 

The exclusion limits for selectrons decaying via a A' coupling are shown in Figure 10(b). 
The lower mass limits for selectrons decaying directly or indirectly via a A or a A coupling are 
included in Table 12. The mass limits for the indirect decays are quoted for a low-mass Xi 
(10 GeV). 
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Figure 7: Selectron indirect decays via a A coupling: upper limits at the 95% C.L. on the pair- 
production cross-sections of a right-handed decaying indirectly. The dashed line shows the 
lowest upper limit (errverTv final state) while the solid line shows the highest one (cfifiuefifii' 
final state). 
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Figure 8: Charged slepton direct decays via a X' coupling: upper limits at the 95% C.L. on the 
pair-production cross-sections of e (solid line) and jl/f (dashed line) decaying directly. 
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m(eR) (GeV) 

Figure 9: Selectron indirect decays via a X' coupling: upper limits at the 95% C.L. on the pair- 
production cross-sections the of a decaying indirectly in the electron channel (solid line), in 
the muon channel (dashed line) and in the tau channel (dotted line). 



Charged Slepton Lower Mass Limits(GeV) 


Species 


A 


A' 




Direct 


Indirect 


Direct 


Indirect 


e 


89 


99 


89 


92 


/i 


74 


94 


75 


87 


f 


74 


92 


75 





Table 12: Lower mass limits for charged sleptons decaying directly or indirectly via a X or a 
X coupling. The mass limits for indirect decays are quoted for a low-mass Xi (10 GeV). The 
limits refer to right-handed charged sleptons for the indirect decays and to left-handed charged 
sleptons for direct decays. 
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Figure 10: CMSSM exclusion region for e+e production in the (me,m^o) plane at 95% C.L. 

for (a) a X coupling and (h) a X' coupling. The limits for X' couplings are presented separately 
according to the first index. For the direct decays, the exclusion is shown for the only possible 
case ofe^eL. The kinematic limit is shown by the dashed lines. 



10.2 Smuon limits 



Upper limits on the cross- sect ions of pair-produced smuons decaying directly via a A coupling 
are shown in Figure 6. 

Figure 11 shows upper limits on the cross-sections of pair-produced smuons followed by an 
indirect decay via a A coupling. 

Upper limits on the cross-sections for pair-produced smuons decaying indirectly via a A' 
coupling are shown in Figure 12. 

Upper limits on the cross- sect ions of pair-produced smuons decaying directly via a A' cou- 
pling to a four-jet final state are shown in Figure 8. 

The exclusion limits for smuons decaying via a A coupling are shown in Figure 13(a). The 
exclusion limits for smuons decaying via a A coupling are shown in Figure 13(b). The lower 
mass limits for smuons decaying directly or indirectly via a A or a A coupling are included in 
Table 12. 



10.3 Stau limits 

Upper limits on the cross-sections of pair-produced staus decaying directly via a A coupling are 
shown in Figure 6. 

Upper limits on the cross-sections of pair-produced staus decaying indirectly via a A coupling 
are shown in Figure 14. 
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Figure 11: Smuon indirect decays via a A coupling: upper limits at the 95% C.L. on the pair- 
production cross-sections of right-handed fiR decaying indirectly. The dashed line shows the 
lowest upper limit (fiTTUfiTTu final state) while the solid line shows the highest one (fififiufififiu 
final state). 
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Figure 12: Smuon indirect decays via a X' coupling: upper limits at the 95% C.L. on the pair- 
production cross-sections of ftn decaying indirectly in the electron channel (solid line), in the 
muon channel (dashed line) and in the tau channel (dotted line). 
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Figure 13: CMSSM exclusion regions for fi~^fi production in the {rrifi, m^o) plane at 95% C.L. 
for (a) a X coupling and (b) a \' coupling. The limits for X' couplings are presented separately 
according to the flavour of the lepton in the decay. For the direct decays, the exclusion region 
is shown for the case flifi-L- The kinematic limit is shown by the dashed lines. 
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Figure 14: Stau indirect decays via a A coupling: upper limits at the 95% C.L. on the pair- 
production cross-sections of right-handed tr decaying indirectly. The dashed line shows the 
lowest upper limit (tttutttu final state), while the solid line shows the highest one (tixiiutihiv 
final state). 
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Figure 15 shows upper limits on the cross-sections for pair-produced staus decaying indi- 
rectly via a A coupling in the electron channel. Since the neutralino can decay into final states 
with either a charged lepton or a neutrino, limits are set by varying this branching ratio for 
the neutrino final state between and 1. For the stau decay in the electron, muon and tau 
channels, the lowest limits arise in all cases for the neutrino branching ratio of 1. This means 
that the limits for the muon and tau channels are identical to the electron channel limit shown. 
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Figure 15: Stau indirect decays via a X' coupling: upper limit on the pair- production cross- 
sections for the indirect decay of a tr in the electron channel. The indirect decay of a tr in the 
muon channel and the indirect decay of a tr in the tau channel yield identical results. 

Upper limits on the cross-sections of pair-produced staus decaying directly via a A' coupling 
to a four-jet final state are shown in Figure 8. 

The exclusion limits in the (mf , m^o^ plane for staus decaying via a A coupling are shown in 
Figure 16(a). The exclusion limits for staus decaying via a A' coupling are shown in Figure 16(b). 
The lower mass limits for staus decaying directly or indirectly via a A or a A coupling are 
included in Table 12. For indirect decays via a A coupling, no exclusion was possible since the 
theoretical cross-sections were smaller than the experimental limits in all cases. 



10.4 Sneutrino limits 

Upper limits on the cross-sections of pair-produced sneutrinos decaying directly via a A coupling 
are shown in Figure 17. 

Upper limits on the cross-sections of pair-produced sneutrinos decaying indirectly via a A 
coupling are shown for Am < m,>/2 in Figure 18. These limits degrade for Am > mi>/2 due to 
the larger fraction of missing energy and lower energy final state charged leptons. 

Upper limits on the cross-sections of pair-produced sneutrinos decaying directly via a A' 
coupling to a four-jet final state are shown in Figure 19. 

Upper limits on the cross-sections of pair-produced sneutrinos decaying indirectly via a A' 
coupling to a four-jet final state with missing energy are shown for Am < mi;/2 in Figure 20. 
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Figure 16: CMSSM exclusion region for t^t production in the {mf,m^o) plane at 95% C.L. 
for (a) a X coupling and (h) a X' coupling. For indirect decays via a X' coupling, no exclusion 
was possible since the excluded experimental cross-section is always larger than the theoretical 
cross-sections. For direct decays, the exclusion region for tlTl is shown. The kinematic limit 
is shown by the dashed lines. 
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Figure 17: Sneutrino direct decays via a X coupling: upper limits at the 95% C.L. on the 
pair-production cross-sections of sneutrinos decaying directly. The dashed line shows the lowest 
upper limit (final states with four taus ), while the solid line shows the highest one (final states 
with four muons). 
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Figure 18: Sneutrino indirect decays via a A coupling: upper limits at the 95% C.L. on the 
pair-production cross-sections of sneutrinos decaying indirectly, for Am < The dashed 

line shows the lowest upper limit (final states with 4 t's and missing energy) while the solid line 
shows the highest one (final states with 4 f^'s and missing energy). 
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Figure 19: Sneutrino direct decays via a \' coupling: upper limits at the 95% C.L. on the 
pair-production cross-sections of (solid line) and Vp^/vr (dashed line) decaying directly. 
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These limits degrade for Am > due to the larger fraction of missing energy and lower 

energy final state jets. 
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Figure 20: Sneutrino indirect decays via a X' coupling: upper limits at the 95% C.L. on the 
pair-production cross-sections of Ue (solid line) and u^/ur (dashed line), for Am < my/2. 

In the CMSSM, the z/g pair-production cross-section is enhanced by the presence of the 
t-channel diagram. The exclusion limits for Vf, decaying via a A coupling are shown in Fig- 
ure 21(a). For a low- mass Xi5 this excluded region degrades due to a substantially smaller 
branching ratio for the decay v —>■ uxi but also due to a smaller detection efficiency for 
Am > mj>/2. The exclusion limits for z>e decaying via a A coupling are shown in Figure 21 (b). 
The lower mass limits for electron sneutrinos decaying directly or indirectly via a A or a A 
coupling are summarised in Table 13. 



Sneutrino Lower Mass Limits(GeV) 


Species 


A 


A' 




Direct 


Indirect 


Direct 


Indirect 


89 
79 


95 
81 


89 
74 


88 



Table 13: Lower mass limits for sneutrinos decaying directly or indirectly via a X or a X' 
coupling. For the indirect decays, lower mass limits are quoted for a large-mass Xi (60 GeV). 

Some combinations of neutralino and sneutrino masses do not correspond to points in the 
CMSSM parameter space for the specific choice of = —200 GeV and tan/5 = 1.5. Regions 
where this is the case are labelled "Theoretically not accessible" in the figures. 

The theoretical cross-sections for and Vt pair-production are smaller than for Vf,. The 
exclusion limits for v^jvr decaying via a A coupling are shown in Figure 22(a). 

The exclusion limits for V^/Vr decaying via a A' coupling are shown in Figure 22 (b). No 
exclusion is possible for the indirect decays. The lower mass limits for muon/tau sneutrinos 
decaying directly via a A or a A are summarised in Table 13. 
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Figure 21: CMSSM exclusion region for z/gZ/g production in the (m,>^, m^o) plane at 95% C.L. 
for (a) a X coupling and (h) a X' coupling, assuming the CMSSM predicted branching ratio for 
the decay z/g J^eXv The kinematic limit is shown by the dashed lines. 
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Figure 22: CMSSM exclusion region for u^u^, v^-v^ production in the {my.m^o) plane at 95% 
C.L. for (a) a X coupling and (b) a X' coupling, assuming the CMSSM predicted branching ratio 
for the decay z> — > uxi. The kinematic limit is shown as the dashed line. 
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Figure 23: Stop direct decays via a A' coupling: upper limits at the 95% C.L. on the production 
cross-section oft in the electron channel (solid line), the muon channel (short-dashed line) 
and in the tau channel (dotted line). Also shown are the maximum (dashed- dotted line) and 
minimum (long-dashed line) cross-sections predicted by the CMSSM, corresponding to mixing 
angles of rad and 0.98 rad (decoupling limit). 

10.5 Stop limits 

A cross-section limit of 0.03 pb was derived for the pair-production of stops decaying directly 
via A i3fc or A 23k (electron and muon channels), in the mass region 45 GeV < mi <100 GeV. 
The excluded cross-sections are shown in Figure 23 as a function of the stop mass for the 
electron, muon and tau channels. 

If one assumes a stop production cross-section as predicted by the CMSSM, masses lower 
than 98 GeV can be excluded for any mixing angle 61. A cross-section limit of 0.07 pb was 
derived for the pair-production of the stops decaying directly via A 33^ (tau channel), in the 
mass region 45 GeV < < 100 GeV. In the tau channel, masses lower than 96 GeV can be 
excluded for any mixing angle 61. More detailed exclusion limits are given in Table 14. The 
value of the mixing angle 61 = 0.98 rad corresponds to the decoupling of the stop from the 
and gives the lowest stop pair-production cross-sections. 

For the light squark decays via A 'couplings, a cross-section limit of approximately 0.1 pb 
was derived for a squark mass up to ~ 70 GeV. This limit degrades to 0.39 pb in the range 
of the W mass as shown in Figure 24. Assuming the same detection efficiencies for the stop 
quarks as for the simulated light squarks, a limit on the stop mass can be established. If the 
stop production cross-section predicted by the CMSSM is assumed, masses lower than 77 GeV 
can be excluded for any mixing angle 61. 



11 Conclusions 



A search was performed for pair-produced sfermions with i?-parity violating decays using the 
data collected by the OPAL detector at centre-of-mass energies of 189—209 GeV, corresponding 
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Stop Lower Mass Limits (GeV) 


Channels 


Oi 




rad 


0.98 rad 


ti e+q 


100 


98 


ti ii+q 


100 


98 


ti T+q 


98 


96 


ti qq 


88 


77 



Table 14: Stop lower mass limits for the two extreme values of the mixing angle in the electron, 
muon and tau channels as well as in the channel. 



OPAL 




m(q) (GeV) 

Figure 24: Squark direct decays via a X " coupling: upper limits at the 95% C.L. on the production 
cross-section of s quarks decaying directly. Also shown are the maximum (dashed- dotted line) and 
minimum ( dashed line ) stop production cross-sections predicted by the CMSSM, corresponding 
to mixing angles of rad and 0. 98 rad. 
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to a total luminosity of approximately 610 pb~^. Direct and indirect i?-parity violating decay 
modes of I and v via the Yukawa A and A couplings were considered, as well as direct i?-parity 
violating decay modes of t and q via A and A . 

No significant excess of signal-like events was observed in the data. Upper limits on the 
pair-production cross-sections for sfcrmions have been computed assuming that only i?-parity 
violating decays occur. These cross-section limits depend only on the masses of the sfermions 
and not on other SUSY parameters. Mass hmits were derived in the framework of the Con- 
strained Minimal Supersymmetric Standard Model whenever the predicted cross-sections were 
sufficiently large. 
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